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H2O maser and a plasma obscuring torus in the radio galaxy NGC 1052 
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We present multi-frequency simultaneous VLB A observations at 15, 22 and 43 GHz towards the nucleus of the nearby 

radio galaxy NGC 1052. These three continuum images reveal a double- sided jet structure, whose relative intensity ratios 

imply that the jet axis is oriented close to the sky plane. The steeply rising spectra at 15^3 GHz at the inner edges of 

the jets strongly suggest that synchrotron emission is absorbed by foreground thermal plasma. We detected H2O maser 

emission in the velocity range of 1550-1850 km s^^, which is redshifted by 50-350 km s^^ with respect to the systemic 

velocity of NGC 1052. The redshifted maser gas appears projected against both sides of the jet, similar to the HI seen in 

absorption. The H2O maser gas is located where the free-free absorption opacity is large. This probably implies that the 

masers in NGC 1052 are associated with a circumnuclear torus or disk as in the nucleus of NGC 4258. Such circumnuclear 
structure can be the sense of accreting onto the central engine. 

© 2009 WILEY-VCH Vcrlag GmbH & Co. KGaA, Wcinheim 



1 Introduction 

The elliptical galaxy NGC 1052 with a LINER is the near- 
est object among radio-loud and broad megamaser sources 
(e.g. Gabel et al. 2000). The redshift of z = 0.0049 (Knapp 
et al. 1978) corresponds to a distance of 20.3 Mpc and the 
angular scale is 98 pc arcsec^^ if we adopt the Hubble con- 
stant Hq = 72 km s~^ Mpc~^ and qo — 0.5. (Spergel et al. 
2003). This galaxy hosts a well defined double-sided radio 
jet of a few pc with P. A. ~ 65°. The jet emanates from the 
nucleus and can be traced out to kilo-pc-scales (Cohen et 
al. 1971; Wrobel 1984; Jones 1984, Kellermann et al. 1998; 
Kameno et al. 2001; Vermeulen et al. 2003; Kadler et al. 
2004b). The lower limit to the jet incUnation of 57° is esti- 
mated by Vermeulen et al. 2003. Radio observations of the 
low-luminosity active galactic nucleus (AGN) of NGC 1052 
with VLBI at multiple frequencies have revealed the pres- 
ence of a dense circumnuclear structure, which obscures the 
very center of this elliptical galaxy (Kellermann et al. 1998; 
Kameno et al. 2001; Vermeulen et al. 2003; Kadler et al. 
2004b). 

This object shows a convex radio spectrum peaked at 
10 GHz, to be classified as a GHz-Peaked Spectrum (GPS) 
source (de Vries et al. 1997; O'Dea 1998). The convex 
spectum is explained as free-free absorption (FFA) due to 
a dense plasma obscuring the nuclear torus (Kameno et al. 
2001, 2003; Kadler et al. 2004b; Vermeulen et al. 2003). X- 
ray spectra also imply a high column density of 10^^ cm^^ 
to 10^'^ cm~^ toward the center, and support the presence of 
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a dense gas torus (Guainazzi & Antonelli 1999; Weaver et 
al. 1999; Kadler etal. 2004a). 

The centerof NGC 1052 harbors a luminous H2O mega- 
maser, which is redshifted by 50-350 km s ^ ^ with respect to 
the systemic velocity of the galaxy (1491 km s^^; de Vau- 
couleurs et al. 1991). The spectral profile typically shows 
a broad velocity width of --100 km s^^ (FWHM) (Braatz, 
Wilson & Henkel, 1994, 1996, Braatz et al. 2003). Past 
VLBI images reveal that H2O maser gas with the velocity 
range of 1585-1685 km s^^ is distributed along the con- 
tinuum ridge of the western jet, 0.05-0.1 pc shifted to the 
west from the gap between the eastern and western jets in 
November 1995. Excitation by shocks into the dense molec- 
ular clump which lies in or around the radio jet, or ampli- 
fication of the radio continuum emission of the jet by fore- 
ground molecular clouds, were suggested by Claussen et al. 
(1998). On the other hand, Kameno et al. (2005) presented 
the circumnuclear torus model to explain the time variabil- 
ity of the H2O maser emission. Relevant to the interpre- 
tation of the H2O maser emission line, several absorption 
lines are also found toward the center of NGC 1052 (H I: 
van Gorkom et al. 1986; OH : Omar et al. 2002; HCO+, 
HCN and CO : Liszt and Lucas 2004). 

2 Observations and Results 

In order to confirm the positional relation between the H2O 
maser gas and the proposed circumnuclear torus, we ob- 
served the continuum and maser emissions in the nucleus 
of NGC 1052 with the VLBA. The continuum emissions 
were measured at 15 and 43 GHz. At 22 GHz, the H2O 
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Fig. 1 The aligned images of NGC 1052 at 15, 22 and 43 
GHz with the VLBA. The synthesized beam sizes (FWHM) 
are 1.3x0.49mas inP.A.= -4.8°, 0.86x0.32 mas inP.A.= 
-7.1°, 0.39x0.17 mas in P.A.= -2.8°, respectively, as 
shown in the left of each image. The contours start at So- 
level, increasing by a factor of 2, where a = 0.24, 1.07 and 
1.45 mJy beam^^, respectively at 15, 22 and 43 GHz. 



maser emission and the line-free continuum emission were 
observed simultaneously. Channel maps of H2O maser 
emission were made every 6.74 km s^^, averaging every 
10 spectral channels. 

The continuum images at 15, 22 and 43 GHz show the 
two-sided jet structure which consists of several compo- 
nents (Fig.[T]l. The extended structures of knots A, CI and D 
become fainter at higher frequency bands. Knots B, C2 and 
C3 are resolved into several components at 22 and 43 GHz. 
The central engine is supposed to lie between the eastern jet 
(knots A and B) and the western jet (knots C3, C2, CI and 
D). The 43 GHz image reveals another knot (C4) located 
between B and C3. The knot C4 is similar to the component 
A15 by Kadler et al. (2004b), which is supposed to be the 
innermost part of the eastern jet. In the 43 GHz image, knot 
A has poorly defined morphology, and knot CI is split into 
several faint peaks. These images at different frequencies in- 
clude uncertainty in absolute positional information through 
the self-calibration process. For alignment of these images, 
we used the relative positions of knots B, CI and C2, which 
are clearly seen in the restored images with a 1.30x0.49 
mas beam at 15, 22 and 43 GHz. Then we derived relative 
offsets to minimize the positional residuals (e.g. Kameno et 
al. 2001; Kameno et al. 2003; Sawada-Satoh et al. 2008). 
Finally, we could overlay the image at 15 GHz and 43 GHz 
with the image at 22 GHz, with positional errors of (±0.05, 
±0.04) mas and (±0.05, ±0.09) mas, in R.A. and DEC re- 
spectively. 




Offset along the jet axis [mas] 



Fig. 2 Free-free absorption opacity along the jet axis (P.A. 
= 65°) which is shown as a solid line in Fig. [3] The fit in the 
inner edge has larger errors, because the continuum spec- 
trum of the inner edge does not show the peak between 15 
and 43 GHz. 

3 Discussions 

3.1 Free-Free Absorption Opacity Distribution 

After restoring with the same beam, we obtained spectral in- 
dices along the jets. The spectral indices indicate that most 
parts of the two-sided jet structure have optically thin syn- 
chrotron spectra at 15-43 GHz except at the inner edge of 
the inner part of the both jets; a steeply rising spectrum 
(a^i = 3.2 ±0.1, a% = 3.1 ±0.1; S cx v<^) at the western 
edge of knot B and at the eastern edge of knot C3 are re- 
vealed. The spectral index exceeds the theoretical limit for 
synchrotron self-absorption (a = 2.5). The highly rising 
spectrum of the inner edge of the both jets implies that the 
synchrotron emission is obscured through the free-free ab- 
sorption (FFA) by the foreground dense plasma, and this is 
consistent with past multi-frequency observations (Kameno 
et al. 2001; Vermeulen et al. 2003; Kadler et al. 2004b). 

Fitting the continuum spectrum at 15^3 GHz to FFA 
model iSu oc i^^" exp (— rz^^^-^)), we obtained FFA opac- 
ity (r) distributions along the jet axis (Fig.|2]l, which reveals 
that high opacity (r > 1000) is found in the inner edge of 
the jets. It implies that the dense cold plasma cover ~2 mas 
in the inner edge, where the central engine is supposed to 
exist. 

3.2 H2O maser emission 

Significant H2O maser emission within the velocity range 
of 1550-1850 km s^^ were detected, 50-350 km s^^ red- 
shifted from the systemic velocity of the galaxy, and it is 
consistent with past single-dish observations (Braatz et al. 
1994, 1996; Braatz et al. 2003; Kameno et al. 2005). The 
maser spots consist of two clusters; the eastern cluster and 
the western cluster are located on knots B and C3, respec- 
tively. The H2O masers projected on the approaching jet, 
or the eastern cluster were detected for the first time. In the 
western cluster, maser spots are distributed along the jet axis 
direction that span ^1 mas (0.1 pc), and show some velocity 
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Fig. 3 Relative distributions of H2O maser spots (full cir- 
cle) with respect to the continuum image at 22 GHz (con- 
tour). The inset box shows the detailed distributions of 
maser spots. 
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Fig. 4 Position- Velocity diagram of H2O maser spots 
along the jet axis with respect to the peak of knot B. A clear 
velocity shift as a function of relative offset along the jet 
axis (P.A.=65° in Fig. |2]) is seen in the western cluster (b). 
On the other hand, the eastern cluster seems to include sev- 
eral sub structures in velocity (a). 



shift along the jet axis direction. Position-velocity diagrams 
of the H2O maser spots along the jet axis (Fig. ^) also re- 
veal the trend of velocity gradient of ^ 250 km s^^ mas^^ 
in the western cluster. We note that Claussen et al. (1998) 
showed a velocity gradient along the east-west direction of 
magnitude ^100 km s^^ mas^^ in the maser cluster on the 
western jet knot. A velocity gradient in the eastern cluster is 
not obvious (Fig.|4^). 

3.3 Interpretation of the nuclear structure 

The brightness temperature ratio between the approaching 
and receding jets is related to the viewing angle of the jet 
axis and the true jet speed. Assuming that knots A and CI 
form a symmetric pair of knots on either side of the nucleus, 
their intensity ratio R is given by 

3-a 

(1) 



R — 77T 



1 + /3 COS 6* 



1 



P cos 9 ^ 

where and Tl^^ are the brightness temperature of knots 
A and CI, respectively, /3 is the true jet velocity as a fraction 
of the speed of light (v /c), 9 is the viewing angle of the jet 
axis, and a is the spectral index. Adopting = 0.64 (Kadler 
et al. 2004b) and a — —1, the viewing angle is estimated to 
be 79-80° and 76-90° using the flux density at 15 and 22 
GHz, respectively. Kadler et al. (2004b) obtained the mini- 
mum value of (3 = 0.64 by taking the jet inclination derived 
by Vermeulen et al. (2003). The jet axis is considered to be 
nearly parallel to the sky plane. 

There are two main ideas about what the H2O masers 
of NGC 1052 are associated with, the jet or the circumnu- 
clear torus (Claussen et al. 1998; Kameno et al. 2005). Here 
we discuss the both ideas to account for the results from the 
observations. One possible scenario is that the H2O masers 
are associated with the circumnuclear torus with the X-ray 
dissociation region (XDR). The H2O masers are distributed 
where the FFA opacity is large, and this suggests that the 
H2O masers and the plasma exist close to each other. In 



the case of the H2O megamaser emission in NGC 4258, 
Neufeld & Maloney (1995) and Herrnstein et al. (1996) pro- 
posed that the molecular disk consists of several layers in- 
cluding a heated molecular layer where the H2O masers re- 
side. 

Kameno et al. (2005) applied this idea to the circumnu- 
clear torus model in NGC 1052. A hot (- 8000 K) plasma 
layer is created on the inner surface of the torus because of 
the direct exposure to the X-ray radiation from the central 
source. This layer is responsible for the free-free absorp- 
tion. The XDR which lies immediately next to the plasma 
layer inside the torus, is heated above ~ 400 K as it is still 
partially irradiated by the X-ray radiation (Maloney 2002). 
Excited H2O molecules in the XDR will amplify the contin- 
uum seed emission from the jet knots in the background and 
result in maser emissions. The presence of masers on both 
jets indicates the thickness of the torus along the orientation 
of the jet, covering at least knot B and C3. If the orientation 
of the jet axis is parallel to the sky plane, the thickness of 
the torus should be therefore 0.2 pc at least. More dominant 
FFA and H2O masers on the receding jet support the cir- 
cumnuclear torus scenario, since the path length within the 
torus toward the receding jet would be greater 

The FFA opacity ris a function of the electron density 
rie and temperature along the line of sight as 



T = 0.46 



2rTi— 1,5 



dL 



(2) 



LOS 



The maximum opacity t ^ 1000 in the path length of ^ 
1 pc gives TigTg^^'^ — 2100. We can estimate the electron 
density Ue > 4.6 x 10^ cm~'^ using Ionization condition 
Te > 10'' K. The electron column density UeL would be 
1.4xl0^^(Te/10''K)°-^^ cm~2, and the value is comparable 
to the atomic column density (~ lO^'^ cm^^) obtained by 
ROSAT and ASCA X-day observations (Guainazzi &An- 
tonelli 1999). 

If the H2O masers are associated with the torus, the red- 
shifted spectrum of the H2O maser emission would be ac- 
counted for by a contraction toward the central engine. The 



www.an-journal.org 



© 2009 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



252 



S. Sawada-Satoh et al.; A Dense Plasma in NGC 1052 




Fig. 5 A cartoon showing the possible environment in the 
circumnuclear torus and jets in NGC 1052. Double-sided jet 
axis inclined by > 76° with respect to the line of sight. In- 
ner surface of the torus is ionized by X-ray emission. XDR 
is formed on the inner layer of the torus and amplifies back- 
ground continuum emission from the jet knot. On B and C3, 
we can see the H2O maser emission and FFA absorption as 
well. On the other hand, only FFA appears on C2. Since 
the gas inside the torus is falling toward the central engine, 
redshifted H2O maser emission is detected. 

positional-velocity diagram along the jet axis for each maser 
cluster (figure m gives the appearance that the H2O maser 
gas closer to the central engine is more redshifted. Such a 
velocity shift as a function of positional offset could indi- 
cate the acceleration of the infalling gas toward the central 
engine. Thus, the circumnuclear torus scenario can explain 
the observed characteristics. 

Alternatively, the H2O masers could be excited by the 
interaction between the jet and circumnuclear molecular gas 
as Claussen et al. (1998) proposed, and as interpreted in 
the case of Mrk 348 (Peck et al. 2003). A weak point of 
the jet excitation scenario is the difficulty to explain why 
the eastern cluster of H2O masers projected against the ap- 
proaching jet is also redshifted. For the masers in the west- 
ern cluster moving with the most redshifted velocity 
400 km s^^) from the systemic velocity of the galaxy, the 
maser gas should have moved 1.9 x 10^"^ pc eastward from 
November 1995 to July 2000, because the jet axis is close 
to the sky plane. This motion is too small to detect for the 
five-year multi-epoch observations. Further VLBI observa- 
tions are necessary in order to detect the proper motion of 
H2O maser spots predicted by the jet excitation model, in 
order to support or refute this scenario. 
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